Enhancement of Andreev reflection by spin polarization in 
out-of-equilibrium ferromagnet-superconductor junctions 
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We report on transport measurements on Nb/Al/Gd/Al/Nb junctions. Bulk Gadolinium is a 
weakly polarized ferromagnet (5-7%), and is present in the junction in granular form (superpara- 
magnet). We show that Andreev reflection is strongly enhanced by a weak polarization, obtained 
by applying an external magnetic field. A new model is proposed that accounts for this effect. The 
transport is described in terms of "hot" carriers that experience Zeeman splitting due to the spin 
polarized background. The Landauer formula is G = H-A(hb), hb being the exchange field at 
the grain boundary, and A(E) the Andreev current coefficient, that increases with E for resistive 
junctions. 
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Much interest has been recently devoted to understand 
the interplay between magnetism and superconductiv- 
ity. For example, recent experiments on ferromagnet- 
superconductor superlattices have shown anomalous os- 
cillating critical temperatures [QJ|] due to a possible cou- 
pling between the magnetic layers through the super- 
conductor. Other experiments have shown new coher- 
ence phenomena in superconducting-ferromagnetic sys- 
tems [||-[| directly related to the presence of magnetic 
ordering. The problem which arises is to clearly un- 
derstand the behavior of normal metal-superconductor 
(N/S) junctions in the presence of ferromagnetism and 
more specifically, the importance of spin polarization on 
Andreev reflection in such junctions. In a N/S sandwich, 
the zero temperature current is transmitted below the 
superconducting gap via Andreev reflection only : an in- 
coming electron with a spin a is reflected at the interface 
as a hole in the band of opposite spin — a, whereas a 
Cooper pair is transferred into the superconductor 
Blonder, Tinkham and Klapwijk (BTK) have shown 
that Andreev reflection is weakened when interfacial elas- 
tic processes are taken into account under the form of a 
repulsive potential V(x) — H5(x) at the interface. This 
theory interpolates between a perfect transparent inter- 
face and an insulating barrier as the dimensionless barrier 
strength Z — H /Tivf increases from zero to infinity. It is 
of importance to understand the behavior of such junc- 
tions when the normal metal is ferromagnetic. In this 
case, assuming thermal equilibrium within the ferromag- 
net, one spin population is depleted with respect to the 
other. Considering channels in the spin up band, and 
Ni < N-\ channels in the spin down band, de Jong and 
Beenakker @] noticed that only N± channels are available 
for Andreev reflection . The corresponding Landauer for- 
mula is G = 4A r |e 2 jh. Andreev reflection is thus strongly 
dependent on the Fermi surface spin polarization. 

Recent experiments in ferromagnet-superconductor 



(FM/S) junctions by Soulen et ai. || (with a ferro- 
magnetic metal FM=Nio. 8 Fe .2, Ni, Co, Fe, NiMnSb, 
Lao.7Sr .3Mn03) and by Upadhyay et al. [jl0| (with 
FM=Co, Ni), have proved that Andreev reflection is 
strongly suppressed as the Fermi surface polarization is 
increased, in full agreement with the predictions by de 
Jong and Beenakker ||. For instance, the low temper- 
ature and low voltage differential conductance Gs(0) = 
(dG(V)/dV)/G n in N/S junction is twice its value above 
the superconducting gap when N is a non magnetic ma- 
terial, because of Andreev reflection. On the other hand, 
Gs(0) falls to zero in the most polarized FM/S junctions 
(90 % polarized if FM=Cr02), showing the suppression 
of Andreev reflection by spin polarization |^|nj . 

The aim of the present Letter is to show that spin 
polarization can enhance Andreev reflection in S/FM/S 
junctions with a sizeable interfacial scattering [Z w 2). 
This provides the opportunity of experimenting FM/S 
junctions in between a point contact regime with Z » 
H|o| and a tunneling regime with Z ^> 1 |riTj . We used 
Gadolinium (Gd) as ferromagnetic metal, in which the 
conduction band is weakly polarized, between 5 and 7 % 
p^^3| ; a small value compared to the polarization of 
37 %in e.g. NiFe reported in the less polarized sample 
studied in Ref. |^|. Unlike the transition metals with s-d 
itinerant magnetism, Gd has 4f localized moments, fer- 
romagnetically aligned through the RKKY interaction 
mediated by the conduction electrons. In spite of the 
complexity of the conduction band structure, we believe 
that the physics involved in our experiment can be cap- 
tured by assuming a Stoner model for the conduction 
band of width W wlO eV @ in the presence of an ex- 
change field /led of the order of a 140 meV for bulk Gd 

& 

A magnetic field applied parallel to the junction is used 
to tune the spin polarization. We show that even with 
polarizations as small as 5 — 7% in the conduction band 
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of bulk Gd, the low voltage resistance of the junction is 
strongly dependent on its polarization. More specifically, 
the resistance of the junction decreases when an external 
magnetic field is applied. We demonstrate that this be- 
havior originates from an enhancement of Andreev reflec- 
tion when the boundary spin polarization of the granular 
Gd increases. We set up a new model which accounts for 
our observations. It is based on the transport of "hot" 
electrons in the Gd grains, which do not thermalize with 
the background ferromagnet conduction band, but expe- 
rience Zeeman splitting. This model relies on a reduction 
of the Gd magnetism at the grain boundaries. 

We study the specific junction Nb/Al/Gd/Al/Nb. The 
layers are prepared by e-beam evaporation in ultra-high- 
vacuum on a silicon substrate held at room temperature. 
Layer thicknesses are controlled during growth with a 
quartz crystal monitor. A typical rate of 1 A/s is used 
for Gd. The base pressure is below 3 x 10 -9 Torr, and 
the working pressure is 2 x 10 -9 Torr for Gd, 4 x 10~ 9 
Torr for Al and Nb. The junction area (0.1 x 0.1 mm 2 ) 
is defined by evaporating insulating SiO through shadow 
masks at 2x 10~ 8 Torr typically. The inset of Fig. [j] shows 
a schematic cross section of the junction, which is made 
of two Nb leads of 250 A, and a sandwich Al/Gd/Al with 
a Gd thickness ranging from 20 and 100 A. A 1500 A Al 
layer has been used to locally lower the superconducting 
T c of the junction and avoid spurious effects due to the 
Nb leads transitions. 
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FIG. 1. Current- voltage characteristics at T=2 K with a 
thickness of 40 A Gd. The top left inset shows the modulation 
of the critical current under low magnetic field. A schematic 
sample cross section is shown on the bottom right inset. All 
samples with different Gd thicknesses exhibit a similar I-V 
characteristic. The magnetic field is applied in the plane of 
the junction. 

Several methods were used to characterize the Gd 
layer. First, Transmission Electron Microscopy images 
on Al/Gd reference samples evaporated under the same 
conditions as the junctions, have shown that Gd lay- 
ers thinner than 20 A are not continuous, thus impos- 



ing a lower bound on samples thicknesses. Gd is poly- 
crystalline and the size of grains ranges between 70 to 
100 A. For the smaller thicknesses, the grains are elliptic, 
of typical dimension 20 Ax 70 Ax 70 A. The magnetic 
properties of Gd were investigated by means of SQUID 
magnetometry. In Fig. ||, we show a Zero Field Cooled 
(ZFC) and Field Cooled (FC) magnetization curve in a 
low applied field (11=50 Oe). These curves exhibit a large 
thermomagnetic hysteresis with an irreversibility temper- 
ature (Tj) of 150 K, and a blocking temperature of about 
120 K. We estimate that the Curie temperature range 
between 50 K and 100 K. This behavior is characteristic 
of a superparamagnet, in agreement with the granularity 
of the magnetic layer. From the saturation magnetiza- 
tion, the magnetic moment per atom is estimated to be 
4/^s, substantially below the Gd bulk moment of 7.6[1b- 
We attribute this reduction of the saturation magneti- 
zation compared to the bulk value to frozen spins that 
do not align in the applied magnetic field due to spin 
glass effects as it is expected for Gd small grains [|l6|,0 • 
The coercitive field of the Gd layers is around 100 mT, 
for all the Gd thickness. The macroscopic magnetization 
of the granular Gd increases as the external field is in- 
creased because the giant spin grains align on average in 
the external field. 
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FIG. 2. Magnetization of a 100 A Gd layer versus temper- 
ature measured in a ZFC-FC procedure in low magnetic field 
H=50 Oe. The Nb layer is not present in this measurement. 
Inset : hysteresis loop of the same sample at 2 K. 

All the transport measurements have been performed 
with a highly sensitive dc-current method using a four 
probe technique |18[ . Resistances range from a few mf2 
to a few /if2 and the voltage across the junction is of 
the order of a nanovolt. A typical current- voltage I-V 
characteristics is shown on Fig. [|. It displays a sizeable 
critical current I c , whose amplitude strongly decreases 
with the Gd thickness, indicating a Josephson coupling 
between the superconducting electrodes through the fer- 
romagnet. This will be exposed in details in a forthcom- 
ing publication. The modulation of I c in a weak applied 
magnetic field (inset of Fig. [j]) confirms the existence of 
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a Josephson coupling and excludes that pinholes domi- 
nate the transport properties. In zero magnetic field (see 
Fig. ||) a first transition appears in transport measure- 
ments corresponding to the superconducting transition of 
Nb between 7.5 and 8 K, that induces superconductivity 
in Al by proximity. At lower temperatures the resistivity 
vanishes due to Josephson effect (inset of Fig. ||). Fig. || 
shows the effect of a magnetic field on the junction re- 
sistance, the magnetic field being applied along the layer 
and perpendicular to the current. The resistance is dras- 
tically reduced by the external magnetic field. We can 
exclude that this behavior originates from part of the bi- 
layer Nb/Al becoming normal. The parallel critical field 
{Hi i) in a thin film is H/ 1 = \/2A^^-H c , where A(<) is the 
penetration depth in the superconductor and d the film 
thickness Q. For Niobium this leads to ff//Nb = 1-7 T, 
a value much larger than the fields used in this investi- 
gation. As for Al, the magnetic field required to destroy 
superconducting proximity effect is larger than the bulk 
critical field H c2 « 200 mT of Nb @. Therefore, the 
bilayer remains superconducting for a magnetic field up 
to 100 mT. 



over, the hysteresis in the magnetoresistance clearly in- 
dicates the onset of a magnetic effect in the transport in 
the junction (inset of Fig. ||). 
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FIG. 3. Ohmic resistance versus temperature measured 
with a bias current larger than twice the critical current under 
a magnetic field applied parallel to the layers (a) H=0 mT, (b) 
40 mT, (c) 60 mT, (d) 100 mT. The solid lines show the fits to 
the N/FM/S model. We have used Z = 2.3 and h b = 0, 0.38, 
0.61, and 0.83 respectively, in units of the Al superconduct- 
ing gap of the order of a meV. The inset shows the behavior 
of the transition temperature under different bias currents in 
zero external field. 



Fig. ^ shows the variations of the ohmic resistance at a 
fixed temperature as a function of the applied magnetic 
field. The sample is first zero field cooled from 10 K to the 
measurements temperatures. The magnetic field is swept 
from -150 mT up to 150 mT, and back to -150 mT. The 
crossover field for the negative magnetoresistance is in- 
dependent on temperature, which is an additional proof 
that this effect does not originate from the destruction 
of superconductivity in the Nb/Al bilayer, in which case 
this crossover field would depend on temperature. More- 
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FIG. 4. Negative magnetoresistance obtained for different 
temperatures on a sample with a 40 A Gd layer and a bias 
current of 5 mA. The details of the hysteresis at 3 K are shown 
in the inset, the arrows indicating the direction in which the 
field is swept. 



The de Jong-Beenakker theory cannot account for the 
present observations, since it would predict an increase 
of resistance upon spin polarization, and also because the 
relative change in resistance induced by a small spin po- 
larization of 5-7 % would be of the order of the spin polar- 
ization pi] ] , incompatible with the magnitude of the resis- 
tance variations shown on Fig. |[ We believe that the fail- 
ure of the de Jong-Beenakker theory in our system lies in 
the fact that transport in the Gd grains is dominantly an 
out-of-cquilibrium transport since the grains (70-100 A) 
and the Gd thicknesses are small. We present here a mod- 
eling of this situation that qualitatively accounts for the 
negative magnetoresistance. In this model, the charge 
carriers are "hot" electrons, which wave vector and spin 
remain constant while they cross the grains. Given the 
small size of the grain, we neglect any inelastic scattering. 
The spin energy of an electron (hole) is — (+)/iGd(x)CT, 
^Gd(x) being the space-dependent exchange field. The 
relevant quantity for Andreev reflection is the exchange 
field /ib at the grain boundary, a quantity which can be 
much smaller than the bulk value had — 140 meV. The 
magnetism of Gd is very sensitive to growth parameters 
p2| , p3| , and surface magnetization can be strongly weak- 
ened by possible adsorbtion of oxygen atoms at the edge 
of Gd grains p4| . An applied magnetic field will favor 
ferromagnetism at the boundaries and increase hb- As 
the external field increases, the nucleation of boundary 
polarization is expected to favor the reversal of the bulk 
giant moment if the grain had a misoriented magnetiza- 
tion. This Gd grain modeling also allows to understand 
that the hysteresis in resistance is weaker than in magne- 
tization, since the former involves boundary magnetism 
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and the latter the superparamagnetic behavior of the gi- 
ant spin grains. 

Now, let us consider the semiclassical treatment of a 
N/FM/S junction ^5|, and assume a scattering only at 
the FM/S interface I2. The presence of a scattering at 
both interfaces leads to similar qualitative conclusions, as 
well as the S/FM/S junction, in which cases a numerical 
treatment is necessary and will be presented elsewhere 
pl| . We note f^ L ^(E, ^1,(2)) the non-equilibrium distri- 
bution functions of right (left) movers with a spin a at 
the Ii=N/FM (I 2 =FM/S) interface. The Zeeman split- 
ting in the FM region is taken into account by proper 
boundary conditions, e. g. f^'(E,Ii) = f T (E + h h a) 
at the Ii interface, fx being the equilibrium distribution 
function in the reservoirs. The exchange field and the 
superconducting gap are assumed to have step function 
variations at the interface I2. We have point split the 
two order parameters at this interface by assuming that 
electrons incoming from the S part will first experience 
the barrier potiential and next Zeeman splitting, or the 
opposite. In the first case, the barrier conductance is 
independent on spin polarization, whereas in the second 
situation, in which both orders coexist close to the inter- 
face, the FV relation is 

/ = \ J (1 + A{E) - B(E)) [f T (E -eV- h h ) (1) 
+f T {E -eV + K) - f T (E - h h ) - f T (E + /i b )] , 

with A{E) and B(E) are the Andreev and backscattered 
transmission coefficients |7j . The Landauer formula is 

G = (4e 2 //i) A(h h ), 

a quantity that increases with the spin polarization hh 
if Z pa 2 0. As shown on Fig. the best fits are in 
a qualitative agreement with experiments. Because of 
the ~ 5 mT hysteresis in conductance, the spontaneous 
boundary magnetization is expected to be finite. How- 
ever we cannot resolve magnetic fields below ~ 5 mT 
whithin the present model. 

Given the weakening of ferromagnetism at the grain 
boundaries |p2| , p3|l , we have shown that the observed 
negative magnetoresistance originates from an enhance- 
ment of Andreev reflection due to Zeemann splitting of 
"hot" carriers. This requires the coexistence of super- 
conductivity and ferromagnetism close to the interface. 
A more complete understanding of this interface would 
involve a treatment of the influence of the ferromagnet 
onto the superconductor, and, in view of the reduction 
of magnetism, the possibility of cryptomagnetic p6| , or 
cryptomagnetic-like structures |27|] at the grain bound- 
aries. 
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